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Abstract. It is known that the fluctuations of suitable linear statistics of Haar distributed 
■ elements of the compact classical groups satisfy a central limit theorem. We show that if the 

O ^ corresponding test functions are sufficiently smooth, a rate of convergence of order almost 

1/n can be obtained using a quantitative multivariate CLT for traces of powers that was 
recently proven using Stein's method of exchangeable pairs. 

1. Introduction 

For n G N let M n denote a random n x n matrix, distributed according to Haar measure on 
one of the compact classical groups, i.e. the unitary, orthogonal, and (if n is even) unitary 
symplectic group. Its eigenvalues A n i, . . . , \ nn lie on the unit circle line T of the complex 
plane C. Write L n = h n (M n ) for the empirical measure ~ YTj=i °f the ^nj- ^ * s wen 
known (see [MM 12— | for a recent quantitative version) that, as n — > oo, L n tends a.s. weakly 



c3 



£> . to the uniform distribution on T. Furthermore, for suitable test functions / : T — > R, the 
fluctuations n(L n (f) — EL n (/)) have been proven to tend to a Gaussian limit as n — > oo (see, 
e.g., |DE0 1 [ ISosOO [ I Wie02] ) . In fact, these results are part of a broader interest in fluctua- 
tions of linear spectral statistics of various random matrix ensembles, both from the physics 
(see the seminal paper [CL95J) and mathematics points of view (see, e.g., |Cha09] and the 
^ | references therein). 



For Haar distributed matrices from the compact classical groups, and in the special case 
that / is a trigonometric polynomial, Johansson has obtained in |Joh97] an exponential rate 
of convergence to the Gaussian limit, using sophisticated analytic tools related to Szego's 
strong limit theorem for Toeplitz determinants. In this note we are concerned with the speed 
of convergence in the case of test functions whose Fourier expansion does not necessarily ter- 
minate. In this more general setting we have been unable to find quantitative convergence 
results in the literature. We point out that for sufficiently smooth test functions a rate close 
to 0(l/n) is an easy consequence of a recent quantitative multivariate CLT for a vector 
of traces of powers of a Haar distributed element of the compact classical groups, where, 
crucially, the length of the vector may grow with the matrix size n. This CLT was proven 
by the present authors in [DS11], using the "exchangeable pairs" version of Stein's method 
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(see |CGS11] for background) and building upon a construction of an exchangeable pair that 
had already been successfully used by Fulman [FullO] in the univariate case. Note that a 
different way of using Stein's method to study the linear eigenvalue statistics of different 
random matrix ensembles was devised by Chatterjee in |Cha09j . 

In what follows we concentrate on the case of random unitary matrices. The orthogonal and 
symplectic cases can be treated along the same lines, leading to the same rates of convergence 
towards a (this time not necessarily centered) Gaussian limit. In Section [2] we recall some 
Fourier analysis as well as the crucial quantitative multivariate CLT for traces of powers. In 
Section [3] we state and prove our main result on linear statistics. 



2. Set-up and background 

Let / G L\T) be real valued. We will view / as a 27r-periodic function on R by tacitly 
identifying f(e lx ) with f(x). For j G Z the j-th Fourier coefficient of / is defined as 



1 



2tt 



fj = fU) = tot I f[r) ' ' J " (lr - 



Note that since / is real valued, f(—j) = f{j)- It is a well known fact about Fourier 
coefficients that smoothness of the function implies quantitive information on the decay 
of the coefficients: 

Proposition 2.1. If f e C k (T), then for allO^ j eZ there holds 

ll/ (fc) lli 

1/0)1 < lfji -. 

Proof. See, e.g., |Kat04l 1.4.4]. □ 

Slightly generalizing this, we will consider functions / G L X (T) with the property that there 
exist k > 1 and C K > such that for all ^ j G Z there holds 



I/O') I < 



Proposition 12.11 then says that for / G C fe (T) (jfe > 2) condition flU holds with k = k, and it 
follows from the results in [GralOl Sec. 3.2.2] that (p]) holds with k = 2 if / has a Lipschitz 
first derivative. Note that since \e^ x \ = 1 for all j and x, our assumption k > 1 implies 
that the Fourier expansion of / will converge normally, hence (by compactness of the torus) 
uniformly to /. So we have that 

oo 

E /o)' ;?r = /(o) + E (/o> iis + /O) ' ' ; ") • 

Let M n be a Haar distributed element of U n and write 

1 n 

L n = L n (M n ) = - V5 A 
n z — ' 

i=i 
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where A n i,...,A nn G T are the eigenvalues (with multiplicities) of M n . Observe that by 
Fubini's theorem and the left and right invariance of Haar measure one has 

EL n (f) = h&[f(\ nl ) + ... + f(\ nn )} 



- [ E[f(t\ nl ) + ... + f{tX nn ))dt 

-E / [f{tX nl ) + ... + f{t\ nn ))dt 
n Jt 

n [ f(t)dt 
Jt 



n 



f(t)dt = f(0). 



In view of this, it follows from the above that the n-scaled fluctuation of L n (f) has a pointwise 
expression 

n(L n (f) ~ E(Ln(/))) = /(Am) + • • • + /(A nn ) - n/(0) 



= £ fU) Tr ( M n) + E fU) Tr(Mfl). (2) 

The following quantitative CLT for vectors of traces of powers of Haar unitaries, proven in 
[DS11], will make it possible to control finite sections of n-dependent length of the expansion 
(|2J): Let M = M n be distributed according to Haar measure on U n . For d G N, r = 1, . . . , d, 
consider the r-dimensional complex random vector 

W := W{d, r, n) := (Tr(M d - r+1 ), Tr(M d - r+2 ), . . . , Tr(M d )). 

Let Z := (Zti-r+i ■ ■ ■ , Zd) denote an r-dimensional complex standard normal random vec- 
tor, i.e., there are iid real random variables X d _ r+1 , . . . , X d , Y d _ r+ i, . . . , Y d with distribution 
N(0, 1/2) such that Zj = Xj + iYj for j = d — r + 1, . . . , d. We take E to denote the diagonal 
matrix diag(d — r + 1, d — r + 2, . . . , d) and write := Y}I 2 Z. 

Recall that the Wasserstein distance for probability distributions fi, v on (C , B(C d )) is 
defined by 

<iyv(/-t, v) := sup 1 1 J hdji — J hdv\ : h : C d — > R and ||/i||Lip < 1 j ; 

II ^|| Lip denoting the minimum Lipschitz constant 

\h(x)-h(y)\ tn , 
Lip := sup AM e o, oo] 

of h. If /j. is the law of X and v is the law of Y, we write c?w(X, F) = dw(/x, ^)- 

We are now in a position to state the crucial quantitative CLT for vectors of traces of 
powers |DSllt Thm. 1.1]. 
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Proposition 2.2. Ifn> 2d, the Wasserstein distance between W and Zy, is 



max { (d -r+i)3/2 > (rf-^) 3/2 v^] 



Z E ) = O I I • ( 3 ) 



d w (W, Z S ) = 0[ — ). (4) 



In particular, for r = d we have 



n 

We will also need the following orthogonality relations for traces of powers: 
Proposition 2.3. Letn,i,j e N. Then 

E(Tr(M;)Tr(M^)) = 

and 



E(Tr(M^)Tr(Mn)) = 5ij(j A n). 

Proof. The first equality follows from the unitary invariance of Haar measure, see also |Sto05t 
Thm. 3.29] For the second equality see jDEOll Thm.2.1(b)]. □ 

3. Rates of convergence for linear statistics 
The following theorem is the main result of this note. 

Theorem 3.1. Let M n be a Haar distributed unitary n x n matrix, and let f E L 1 (T) be 
real valued and satisfy condition ([T]) above with k > 1. Then, as n — >■ oo, the fluctuation 
n(L n (f) — KL n (f)) converges in distribution to a centered Gaussian random variable with 
variance 

oo oo 

o 1 : £ |/,| 2 | J | = 2^|/ J | 2 J <oo. 

i=-oo j=i 

If C n denotes the law ofn(L n (f) — EL n (/)), then there exists C = C(f,n) > such that for 
any n E N the Wasserstein distance between C n and N(0, a 2 ) can be bounded as follows: 



fe-i_- 

2 , 



d w (£„,N(0,O) < C n~y^7 

In particular, if f is in C°°(T) ; then dy^{C n , N(0, a 2 )) is O y—r=j) f or a ^ e > 0. 

Remark 3.2. In [DE01J the corresponding weak convergence result without a bound on 
the speed of convergence is proven under the weaker assumption that a 2 be finite, which 
amounts to requiring that / be an element of the Sobolev space if 1//2 (T). Since on the one 
hand, elements of i7 1//2 (T) need not be continuous |Kat04l 1.8.11], and on the other hand, 
our method of proof requires the assumption that k > 1, which implies that / must be the 
uniform limit of a sequence of continuous functions, the present result applies to a strictly 
smaller class of test functions. 
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Proof. It follows from ([T]) that a 2 < oo. Set 



S n := E fU) T rM) + E fO) Tr ( M ^ 

3=1 3=1 

and write, for d = d(n) < | to be chosen later, 

d d 



and 



3=1 3=1 
C(2) ._ C _ q(l) 



(5) 



Then, by Proposition 12.11 and dominated convergence, 



E|^ 2) | 2 = E E (/< Tr «) + /' Tr ( M -)J (/i Tr ( M «) + A Tr ( M ") 



E E (Tr(M;)Tr(M^)J + fj 3 E (Tr(M*) Tr(M^)) 



+ fjj E (Tr(M») Tr(M^)) + fj 3 E (Tr(M* ) Tr(M, 

oo 

= 2 £ |/,f(jAn), 

where we have used Proposition 12.31 in the last step. Consequently, Proposition 12.11 implies 
that 

oo oo 1 

E|^)| 2 <2^i|/ :? .| 2 <2||/W||^ -Lj 
j=d+i j=d+i 3 

POO 

<2||/«||?/ - 

For z G C d define 

d 

<P(?) = Vn{z) = E h Z 3 + h~3- 
3=1 

We write f — (fi, ■ ■ ■ , fd)- Then, for any z,w G C d , Cauchy-Schwarz implies 

d 

E ~ ^) + - ^) 



_ oil fW|| 2 1 

IIJ lll (2A;-2)rf( 2fc - 2 )" 



(6) 



1^0) -^H| = 

We have thus obtained that 



3=1 



<2||/|| 2 ||z-u>|| 2 . 



MlLip < 2||/|| 2 . 



Note that, by Proposition 12. 1[ 



d(n) oo 

£iAi 2 <ii/ w ii?£i<~ 



i=i 



i=i 
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is bounded by a constant independent of n and of the specific form of /. 

Now let C be a real standard normal random variable, and Z, £, Zs as above in Section [2j 
Let g : R — > R be Lipschitz with ||g||Lip < 1- Then 

\E(g(S n ))-E(g{*())\ 

< \E(g(S n ))-E(g(S^))\ + \E(g(S^)) - E{g( v {Z*)))\ + \E{g{ V {Z^))) -E{g{aQ)))\ (7) 
In view of (jSJ) and (jH]), the first summand in (j7j) can be bounded by 



Skip E|Sf | < |M| Lip \/E|s£ 2) | 2 < |M| Lip \\f k \ 



(2k -2) d 2fc " 2 ' 



Writing W n := (Tr(M n ), . . . , Tr(M^)), by Proposition 12.21 the second summand in (J7|) has an 
upper bound 



|E(^W)) -E(^(Z E )))| = \E(g&(W n ))) -E{g{^)))\ 

= \\<p\\u P 



VIMkip / VIMIup 

<IMIii P d w (W n ,Z E )<2||/|| 2 C , 

— n 

where C is the constant from Proposition 12.21 

To bound the third summand in ([7]), we will first study the distribution of <p(Zx). Recall 
that Zy, is the random vector [Z\, y/2Z 2 , . . . , \/dZ d ), where the Zj (J = 1, . . . , d) are iid stan- 
dard complex normal random variables, i.e., there are iid real random variables X\, . . . , Xj, 
Y"i, . . . ,Yd with distribution N(0, 1/2) such that Zj = Xj + %Y 3 - for j = 1, . . . , d. Then 



3=1 j=l 
d 

3=1 
d 

(Re(/ i )X j - ln,(/ ; )V ; 



3=1 

(p(Z-z) is thus a centered real normal random variable with variance 

d . d(n) 

4E|( Re (/i) 2 + M/i) 2 ) = 2^j|/ i | 2 =: a*. 
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So we obtain 



\E(g(<p(Zv))) - E(g(<j())\ = \E(g(cx n O) ~ E( fl K))| 

< \\g\\u P E|a n C - o"C| 

< \\g\\u P Wn - o-| E|c| 

< \\g\\u P Wn - a\- 



Now, 



CO oo 1 

k-^i=2^mi 2 <2|i/ w ii; Ei 

i=d+i j=d+i J 

<nf {k) \\i jf^dx 



211/ 



(fe)|,2 



11 (2Jfe-2) rf 2fc " 2 ' 
Since |cr^ — cr 2 | = |(cr n — cr)(<7„ + a)\, we have that 

|<7 n — (7 1 = < 

u n + o o 
and the third term in (J7|) may thus be bounded by 

?(fc)||2 ' 



Comparing the bounds that we have obtained for the individual summands in (J7J), we obtain 
that 

f 1 d 7/2 1 

|E^(5))-E(^((7C))|<|b||i4pC7 /ffc ™ax| — , — j 



n 



<\\g\\u P c f , k (_ + _), (9) 



where C/^ depends only on / and k. 



One verifies that, on ]0, oo[, the function f(x) := -jkr + 2— has a unique global minimum 
at 

x : = ( J ns+2fe. 

Furthermore, / is strictly falling on ]0, Xq[ and strictly growing on ]x , oof. Setting d\ := \x§\ 
and d,2 '■— d\ + 1, we see that the minimizer of / on the positive integers is in {d\, d 2 }. 

There are < a < (3, depending only on k, such that for all n 6 N one has 

2 2 
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Hence, 



f{x ) < 



i . 4 /2 



< 



i 



d 



a 





a 



where 7 : 



1 



+ /3 7//2 . Comparing with ([9]) yields the theorem. 



□ 
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